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We present a search for the decays B — » K*vV using 454 xlO 6 BB pairs collected at the T(4S) 
resonance with the BABAR detector at the SLAC PEP-II B-Factory. We first select an event sample 
where one B is reconstructed in a semileptonic or hadronic mode with one charmed meson. The 
remaining particles in the event are then examined to search for a B — > K*vv decay. The charged 
K* is reconstructed as K* + -» K^tt + or K* + -> K + ir°; the neutral K* is identified in K*° -» 
K + tt~ mode. We establish upper limits at 90% confidence level of B(B + K* + vV) < 8 x f0~ 5 , 
B(B° -> K*°vV) <12x 10" 5 , and B(B -> K*vV) < 8 x f0~ 5 . 
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PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



In the Standard Model (SM) the b — ► svv process oc- 
curs via one-loop box or electroweak penguin diagrams, 
as shown in Fig. [TJ and it is therefore expected to be 
highly suppressed. Due to the absence of photon pen- 
guin contributions and long distance effects, the corre- 
sponding rate is predicted in the SM with smaller the- 
oretical uncertainties than b — > s£ + £~. In particular, 
the SM branching fraction for B — ► K*vV is expected 
to be (1.3±g;|) x 1CT 5 However, this could be en- 
hanced in many new physics (NP) scenarios, where sev- 
eral mechanisms contribute to the rate. In Ref. [l[ non- 
standard Z° coupling contributions are computed, giv- 
ing an enhancement of up to a factor 10. Moreover, new 
sources of missing energy, such as light dark matter can- 
didates [IJ or unparticles 0, 0] , could contribute to the 
rate and produce a final state with a K* [|| plus miss- 
ing energy. The kinematics of the decay is described in 
terms of s vv = wri„ v /rn B , where m vv is the invariant mass 
of the neutrino pair and m B is the B meson mass. NP 
effects can strongly affect the shape of the s uv distribu- 
tion [l|, |4J , and this is taken into account in the present 
work to obtain a model independent limit. 




b t,c s 




FIG. 1: SM diagrams for b — > svv transitions. 
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A previous search by the Belle Collaboration [5| sets 
upper limits of B(B + — ► K* + vv) < 1.4 x 10~ 4 and 
B(B° -> K*°vv) < 3.4 x 1(T 4 at 90% confidence level [|. 

In this paper we present the first BABAR search for both 
neutral and charged B — > K*vv decays. The analysis is 
based on the data collected with the BABAR detector Q 
at the PEP-II storage ring. The sample corresponds to 
an integrated luminosity of 413 fb _1 at the T(45) reso- 
nance, consisting of about 454 x 10 6 BB pairs. An addi- 
tional sample of 41 fb _1 was collected at a center of mass 
energy 40 MeV below the T(45) resonance in order to 
study continuum events: e + e~ — > qq (q = u,d,s,c) and 
e + e~ — > t + t~. Charged-track reconstruction is provided 
by a silicon vertex detector and a drift chamber operat- 
ing in a 1.5 T magnetic field. Particle identification is 
based on the energy loss in the tracking system and the 
Cherenkov angle in an internally reflecting ring-imaging 
Cherenkov detector. Photon detection is provided by 
a CsI(Tl) electromagnetic calorimeter (EMC). Finally, 
muons are identified by the instrumented magnetic-flux 
return. 

Pairs of photons with invariant mass between 115 and 
150 MeV/c 2 are considered as ir° candidates. The K® 
candidates are reconstructed from pairs of oppositely 
charged pions. 

A GEANT4-based [| Monte Carlo (MC) simulation is 
used to model the detector response and test the analysis 
technique. Approximately 13 x 10 6 events are simulated 
where one B meson decays to a signal candidate mode 
and the other B decay is unconstrained (signal MC sam- 
ple), and the kinematics of the signal decay is described 
by a pure phase space model. Simulated generic BB 
and continuum samples are used to investigate the back- 
ground contamination and perform systematic studies. 

Due to the presence of two undetected neutrinos in the 
final state, the B — > K*vv decays cannot be fully recon- 
structed. Hence, one of the two B mesons produced in 
the T(45) decay (the tagging B) is reconstructed in a 
scmilcptonic (B s \) or a hadronic (Shad) mode containing 
a charmed meson. Then a K* and missing energy are 
searched for in the rest of the event (ROE), defined as 
the set of tracks and EMC clusters not associated with 
the tagging B. The two tagging strategies provide non 
overlapping samples and the corresponding results can be 
combined as independent measurements. Selection crite- 
ria are applied to suppress the background contamination 
and an extended maximum likelihood fit is performed to 
extract the signal yields (N s ), which are finally used to 
determine the decay branching fractions (B). In general, 
these can be written as: 



where e is the total signal efficiency measured with the 
signal MC sample and N BB the number of produced 
BB pairs. In the semileptonic (SL) tagged analysis we 
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adopt Eq. JT]) and use control samples to correct for small 
data/MC disagreements in the efficiency. In the hadronic 
(HAD) analysis, in order to avoid large systematic un- 
certainties associated with the MC estimate of the re- 
construction efficiency for the Shad, we normalize the 
branching fraction with respect to the number of data 
events with a correctly reconstructed -Bhad (Nb 



had ' 



B = 



-BB 



-B 



had 



SB ■ ■ N B , , 
sig "had 



(2) 



- B 



had 



where £s sig is the efficiency related to the signal side re- 



construction and selection, while e 



BB 
3 had 



and e 



K'vV 
3 had 



are the 



Shad reconstruction efficiencies in events with generic BB 
decays and events containing the signal process, respec- 
tively; to account for differences among them, observed 
in the MC samples, their ratio e K R vv le?> B is used in 

"had "had 

Eq. [3 as a correction factor. 

The event selection starts from the reconstruction of 
the tagging B. In the SL analysis, we search for a B — > 
D^lu decay. Neutral D mesons are reconstructed in 
the K~tt + , K~tt + it , K~n + ir + ir~ and K®Tr + Tr~ modes. 
Charged D mesons are reconstructed in the K^ir + Tr + 
and K®ir + final states. The D*° candidates are recon- 
structed in the D*° — > D°j channel and the D* + candi- 
dates in the D*+ -> D tt + or D*+ -> D+tt° channels. Fi- 
nally, a lepton (electron or muon) candidate is associated 
to the D meson and a kinematical fit is performed to find 
the B s \ decay vertex. Preliminary selection requirements 
are applied on the D mass (within 0.07 GeV/c 2 of the 
nominal mass in the K~tt + tt° mode, within 0.04 GeV/c 2 
elsewhere) and the momentum of the lepton in the cen- 
ter of mass (CM) frame (|pf| > 0.8 GeV/c). We also 
require the CM angle between the B s \ and the D^*'l pair 
to satisfy —5.0 < cosOb^di < 1-5, where cos8b,di can be 
calculated from the D^*'l four-momentum assuming that 
only one massless particle is missing: 



cos t) B .Dl 



9 fp* TP* 
^^B^xp^Dl 



m. 



Dl 



2|Pb. 



exp 



IPdjI 



(3) 



In Eq.[3l m B is the nominal B mass, E* Bexp and \Q B ,exp\ 
are the expected B energy and momentum, fixed by the 
energies of the beams and evaluated in the CM frame, 
and |p^>;| is the D^l pair momentum in the CM frame. 
Values of cos 9b, di out of the physical range [-1,1] are due 
to resolution effects and missing particles in the Dl recon- 
struction. The distributions of the D mass and the lepton 
momentum in the CM frame, after the reconstruction of 
the signal B, are shown in Fig. [2l The plots are made af- 
ter the signal reconstruction since in case of multiple B s \ 
candidates, the selection of the best one depends on the 
signal side reconstruction too, as will be discussed later; 
events where in the signal side a K* + — » K + tt° channel 
is reconstructed are shown. If one B — > D^lu candidate 
can be reconstructed in the ROE with the same proce- 
dure adopted for the tag side, the event is selected as 



a control sample of double-tagged events for systematic 
studies. 




1.85 1.9 
m D (GeV/c 2 ) 



1.5 2 
p* (GeV/c) 



FIG. 2: The D mass a) and the CM momentum of the B s \ 
lepton b) , in the SL analysis from MC simulations (histogram, 
the hatched area shows the continuum contribution), on-peak 
data (•) and luminosity scaled off-peak data (o). Events 
where in the signal side a K* + — > K + ir° channel is recon- 
structed are shown. 

In the HAD analysis, we reconstruct Shad decays of the 
type B — ► DY, where D refers to a charm meson, and 
Y represents a collection of hadrons with a total charge 
of ±1, composed of Hitt^ + 112K + n 3 K® + n^n , where 
n x + n 2 < 5, n 3 < 2, and n±_ < 2. Using D a (D + ) and 
D*°(D* + ) as seeds for B~(B a ) decays, we reconstruct 
about 1000 different decay chains. Charmed mesons 
are reconstructed in the same final states used in the 
SL analysis, along with the additional channels D + — > 
K+n+n-n , K^tt+tt+tt-, K q s tt+tt° and D*° D Q ir°. 
-Bhad candidates are selected by the two kinematical vari- 
ables: 



AE 



= K 



7*2 



IpbI 



^bcam 5 



(4) 



where E^ c&m is the beam energy and E B and are 
the energy and the momentum of the Bhad in the CM 
frame. For correctly tagged B candidates, the toes dis- 
tribution peaks at the nominal B mass value and AE at 
zero. Hence, a selection is applied by requiring —0.09 < 
AE < 0.05 GeV and 5.270 < m ES < 5.288 GeV/c 2 . 
The number of correctly reconstructed -Bhad events, to 
be used in Eq. @, is extracted from the toes distri- 
bution of on-peak data. Background events are classi- 
fied in four categories: combinatorial B°B°, combinato- 
rial B + B~ , e + e~ — * cc and e + e~ — * qq (q = u,d,s). 
Other sources of background are found to be negligible. 
For each category, we extract the Toes shape from MC 
simulations. The normalizations of the continuum con- 
tributions are taken from off-resonance data, scaled by 
the luminosity. The normalization of the BB contribu- 
tion is extracted from a \ 2 fit in the 5.22 < toes < 
5.26 GeV/c 2 region. The number of misreconstructed 
Bhad in the signal region is extrapolated from the fit and 
subtracted from the data yield. In Fig. [3] the toes dis- 
tributions for charged and neutral Bhad are shown: the 
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on-peak data are superimposed to the estimated back- 
ground contribution. After background subtraction, in- 
cluding correction factors and systematic uncertainties 
that will be discussed later, we determine NB had = 
(7.175 ±0.008(stat)±0.222(syst)) x 10 5 for neutral B had 



and Nb- 



had 



(10.128 ± O.OlO(stat) ± 0.344(syst)) x 10 5 



for the charged i?had- 




(GeV/c 2 ) 



FIG. 3: The rriES distributions for charged a) and neutral b) 
Shad- The points represent the on-peak data and the hatched 
area shows the estimated background contribution. 

For each reconstructed tagging B, we search for a K* 
candidate in the ROE. A neutral K* can be reconstructed 
in the K + ir~ mode, while a charged K* can be recon- 
structed in the K®ir + and K + ir° channels. The number 
of tracks in the ROE is required to match exactly the 
number of expected tracks for the selected mode. The 
signal B must have opposite flavor (inferred from the K* 
flavor) with respect to the tagging B. 

If more than one -Bsi(-Bhad) — B s i g pair has been recon- 
structed, only one of them is selected. In the SL analysis, 
we adopt a Bayesian approach to define the probability 
that both signal and tag side have been correctly recon- 
structed, given a set x of observed quantities: 



P(TT\x) 



P{x\TT)P(TT) 



; = TT, TF, FT, FF 
(5) 

where TT (FF) indicates that both sides are correctly 
(wrongly) reconstructed and TF (FT) that only the tag 
(signal) side is correctly reconstructed. The candidate 
with the highest P(TT|x) is retained. The set x is com- 
posed of the x 2 probabilities of the B s \ and the K* vertex 
fit. The corresponding likelihoods and prior probabilities 
are modeled from MC simulations with truth information 
to identify the correctly reconstructed candidates. In the 
HAD analysis, if more than one Bhad is reconstructed, 
the best one is selected according to the smallest AE; if 
there are multiple K* candidates associated to the best 
-Bhad, the one with a reconstructed mass closest to the 
world average value [9( is chosen. 

Background contamination is reduced by applying a 
further selection on the BB candidates. Event shape 
variables, namely cos9* BT (the angle between the tag 
side reconstructed momentum and the thrust axis 1 101 ] 
of the ROE) and i?2 (the ratio of the second and ze- 
roth Fox- Wolfram moments [HI]), are used to reject the 



TABLE I: Discriminating variables used in SL and HAD anal- 
yses and specific selection requirements. Values given in the 
squared brackets represent the lower and upper selection cri- 
teria imposed on the respective quantity. 



Variable 


Mode 


Ran 


ge 








SL 


HAD 


pAC f) * _ 


K* + 




[-0 Q8 fl Q71 




K*+ 


-s. K°ir + 


[-0.99,1.00] 


-0.95,0.95] 




K*° 


-> K + TT~ 


[-1.00,1.00] 




Ri 


K*+ 




[0.01,0.82] 






K*+ 




[0.01,0.71] 


0.00,0.70] 




A 


— > l\ TT 


[0.00,0.80] 




TTlK* 


K* + 




fO 83 Q71 




(GeV/c 2 ) 


K*+ 

7 \ 


— > K°TT+ 

— ► J\ S TT 
— > A 7T 


[0.85,0.95] 
[0.84,0.97] 


0.84,0.96] 




K *+ 








(GeV/c 2 ) 


K 
K 


zvO + 

A S 7T T 

TS+ — 

— > K ix 


[0.49,0.50] 


0.49,0.51] 


miss ~r l^miss 

(GeV) 


K *+ 

K*+ 
K*° 


-> K+tv- 


[5.01,8.73] 
[5.11,9.01] 


- 


cos 


K*+ 




[-0.90,0.88] 






K*+ 
K*° 


-> K+n- 


[-0.88,0.85] 
[-0.95,0.89] 


-0.90,0.90] 


'extra 


K*+ 








( GeV) 


K*+ 
K*° 


-^K° s ir+ 

K + TY- 


[0.00,1.20] 




IPfl 


K* + 




[0.95,2.40] 




( GeV/c) 


K'+ 
K*° 


-» K+TT- 


[0.80,2.40] 
[0.84,2.48] 





continuum background. The K* mass (m#-*) and, for 
the K®tt + mode, the K® mass (m K o ) allow rejection of 
combinatorial K* candidates. We define the missing 4- 
momentum due to unreconstructed neutrinos as the dif- 
ference between the T(AS) 4-momentum and the recon- 
structed tagging B and K* 4-momenta. It is exploited 
in the selection through the combination E* niss + Ip^J 
(the sum of the missing energy and the missing momen- 
tum evaluated in the CM frame) and the angle cos 6** niss 
(the azimuthal angle of the missing momentum in the CM 
frame). The extra neutral energy -Eextra, defined as the 
sum of the energies of the EMC neutral clusters not used 



TABLE II: Further selection requirements applied to the B s \ 
candidate (m^ DG is the nominal D mass [s|). 



Variable 


Mode 


Range 


cos 0B,DI 


D° modes 
D modes 


[-2.00,1.00] 
[-1.00,1.00] 


mo - mo 00 
( GeV/c 2 ) 


B+ -> D°(K-n + n°)e + i> 
other B — ► Dlu modes 


[-0.035,0.035] 
[-0.020,0.020] 


Am 

(GeV/c 2 ) 


D*° -» D°j 
D *± _ D±(o) 7r o(±) 


[0.10,0.15] 
[0.14,0.15] 
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FIG. 4: Fit results for the extra EMC energy S e xtra a ) - c ) f° r the SL analysis and the neural network output NN ou t d) - f) 
for the HAD analysis. From left to right, K*+ -» if+Tr , -» K°ir + and A"* -> if + 7r~. Data are shown as points, and 

the fit result is shown with a solid line. The dotted and dashed lines show the estimated signal and background contributions 
respectively. 



to reconstruct either the tag or the signal B, is exploited, 
considering that signal events have no additional neutral 
particles produced in association with the K* . The re- 
quirements applied on the selection variables described 
above are listed in Tab. [TJ 

In the SL analysis, the selection is optimized in the 
MC samples by maximizing the Punzi figure of merit [i~2| , 
given by e/ (n a /2 + ^/N b ), where e is the total signal ef- 
ficiency, N b is the number of expected background data 
events and n a — 1.285 corresponds to a one-side 90% 
confidence level. We also refine the B s \ selection with 
respect to the one applied before the choice of the best 



candidate, and the corresponding requirements are sum- 
marized in Tab. |IIJ where, Am is the difference between 
the D and D* masses, expected to be 142.17 ± 0.07 0. 
The total signal efficiency, evaluated with MC simula- 
tions, is given in Tab. IIIII The variable -Eextra is not 
used in the selection optimization, and its distribution is 
used in an extended maximum likelihood fit in order to 
extract the signal yield. Due to the lower bound on the 
energy of detected photons (50 MeV), the distribution of 
^oxtra is not continuous, so we define the likelihood in 
the following form: 



C(N s ,N b ) 



-[(l-f s )N B +(l-f b )N b ] 



Nil 



[P s ig{E ex tra,i\Vsig){l ~ fs)N s + Pbkg(E e xtra,i\Pbk g )(l - fb)N b ] 

i=l 

(f s N a + f b N b ) N °e-V- N - + f> N >) 
N \ 



(6) 



where N s and N b are the expectation values for the numbers of signal and background events; f s and f b 
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TABLE III: Expected signal and background yields (N s and Nb respectively) from MC studies (assuming the SM B for the 
signal) and results of the data fit, along with signal efficiencies, corrected for systematic effects. Expected signal yields are 
evaluated according to the SM expected B. The first error on the fitted signal yield and on Nb^ , is statistical, the second is 
systematic. The corresponding upper limits are also quoted. 



K* mode 


K+7V [> 




K+tv- 


SL ANALYSIS 


Expected Yields 








N s 


3.31 


2.54 


4.07 


N b 


697 


827 


468 


Scxtra Fit Results 








N s 


-22 ± 16 ± 14 


3 ± 17 ± 15 


35 ± 13 ± 9 


Nb 


754 ± 32 


869 ± 34 


476 ± 25 


e (xHT 4 ) 


5.6 ± 0.7 


4.3 ± 0.6 


6.9 ± 0.8 


N BS (xlO 6 ) 




454 ±5 




UL (90% CL) 


9 x 10" 


■b 


18 x 10~ s 


HAD ANALYSIS 


Expected Yields 








N s 


0.87 


0.77 


1.64 


N b 


46 


35 


73 


NN Fit Results 








N s 


5 ± 6 ± 4 


3 ± 7 ± 4 


-10 ± 9 ± 6 


N b 


39 ± 9 


51 ± 10 


77 ± 13 


e Bsig (XlO"*) 
Nb^„ (xlO 5 ) 


5.8 ±0.5 5.2 ±0.6 
10.128 ±0.010 ±0.344 


16.6 ± 1.4 
7.175 ±0.008 ±0.222 


UL (90% CL) 


21 x 10 


-5 


11 x 10~ 5 



are the fractions of signal and background events with 
-Eextra = 0, and are fixed from the results obtained in 
the MC samples; No and N\ are the numbers of observed 
events with ©extra = and ©extra > 50 MeV respectively; 
and P S ig and Pbkg are the probability distribution func- 
tions (PDF) for signal and background, depending on a 
set of parameters p s j s and Pbkg respectively. MC studies 
show that the background distribution is well described 
by a first-order polynomial PDF, while the signal shape 
can be parameterized with an exponential function and, 
in the charged modes, with an additional Landau con- 
tribution that accounts for photons from a tag side D* 
not associated to the B s \ during the reconstruction. The 
parameters of the PDFs are evaluated in the MC samples 
and fixed when fitting the real data. The fit strategy is 
validated by means of simulation studies which do not 
show any significant bias on the signal yields. The fits to 
the .©extra distributions in the data sample are shown in 
Fig. 2] and the fitted yields are quoted in Tab. IIIII along 
with the total efficiencies e. 

In the HAD analysis, we apply a loose selection 
(Tab. HI, then all discriminating variables are used as 
inputs for a Neural Network (NN), whose output vari- 
able NN ou t is fitted in the region NN ou t > 0.6, where 
the events from the signal MC sample are mostly con- 
centrated. The upper bound of the fit region is differ- 
ent among the three K* modes, reflecting the shape of 
NN out in the signal MC sample. Three different NN 
are trained, one for each K* decay mode. The signal 
output is described with an exponential function for the 



K *o x+TT- mode and a Crystal Ball PDF [3 for the 
charged K* channels. The background is parameterized 
by 

f(x) oc X + k h 1 . (7) 

Also in this case, in the fit to real data the signal and 
background PDF parameters are fixed to the values ex- 
tracted from the MC simulations. Simulated experiments 
are used to validate the fit strategy. The fits to the NN ou t 
distributions in the data sample are shown in Fig. 21 the 
fitted yields and the -B s ; g efficiencies sb ■ are quoted in 
Tab. Hill 

The branching fraction measurement is affected by sys- 
tematic uncertainties related to the signal efficiency esti- 
mate, the B normalization and the signal yield extraction 
from the fit. 

The signal efficiency has an uncertainty due to the 
limited MC statistics. Control samples are used to es- 
timate and correct for data/MC disagreement in the 
charged particle tracking, neutral particle reconstruction 
and particle identification. Uncertainties associated with 
the event selection criteria are computed depending on 
the specific selection strategy. Data/MC comparisons 
and expected detector resolutions provide an estimate of 
possible discrepancies in the distribution of the selection 
variables. For the SL analysis these values are used to 
vary the selection requirements and evaluate the impact 
on the efficiency; for the HAD analysis they are used to 
randomly smear the distributions of the NN inputs and 
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TABLE IV: Summary of systematic uncertainties on the signal efficiency, signal yield, and normalization. 







SL ANALYSIS 






HAD ANALYSIS 




K* mode 


K+n° 


K° s n+ 


K+TY- 


K+ir 


K° S 7V+ 


K+TT- 


Signal efficiency (%) 


MC statistics 


1.4 


1.7 


1.3 


2.9 


3.1 


2.4 


Best pair selection 


0.2 


0.0 


0.0 








Tagging Efficiency 


10.0 


10.0 


10.0 








Tracking 


0.3 


1.0 


0.7 


0.3 


1.0 


0.7 


7T° reconstruction 


3.0 






3.0 






K% reconstruction 




2.5 






2.5 




Particle ID 


1.7 




1.4 


1.7 




1.4 


Selection variables 


5.0 


7.3 


5.1 


5.3 


8.6 


3.8 


Model dependence 


4.5 


4.8 


1.3 


6.3 


7.4 


6.9 


Signal yield (events) 


Signal PDF param. 


0.7 


1.4 


0.2 


0.2 


0.3 


0.2 


Bkgd PDF param. 


11.0 


11.0 


7.7 


2.8 


2.8 


4.5 


Signal PDF shape 








1.2 


1.7 


1.2 


Bkgd PDF shape 


6.4 


4.9 


2.8 


2.1 


1.6 


3.4 


Normalization factor (%) 


N BB or N BhaA 


1.1 


1.1 


1.1 


3.4 


3.4 


3.1 



evaluate the impact on the efficiency after the NN cut. 
The uncertainty due to the residual model dependence 
of our measurement is estimated as follows. We apply a 
weight to each MC event, based on the generated value of 
s„„, in such a way that the weighted distribution matches 
the expected distribution in the SM or some specific NP 
model. Then the efficiency is evaluated taking into ac- 
count the weights and is compared to the nominal effi- 
ciency (obtained from the unweighted MC events, gener- 
ated with a pure phase space model). For the SL analy- 
sis two further uncertainties are associated with the best 
candidate selection and the B s \ selection efficiency. The 
former is evaluated by modifying the input likelihoods of 
Eq. ([5]) according to data/MC comparisons. Concerning 
the B s \ selection efficiency, we apply a correction given 
by the square root of the data/MC ratio of the number of 
double-tagged events. Alternative approaches are used to 
compute the same correction factor and the largest dis- 
crepancy with respect to the nominal approach is taken 
as a systematic uncertainty. 

The error on the number of produced BB events is es- 
timated as described in Ref [lj]. The systematic error 
for Ns had , used in the HAD analysis, is computed by 

varying the MC BB component both in shape and nor- 
malization. The ratio vv Ie?> b , is used to correct 

_ ^had ' ^had ' 

the tag yield and to assign further systematic uncertain- 
ties. Since the tagging efficiency depends on the global 
event multiplicity, this ratio is expected to be different 
from 1 and to depend on the signal side decay and the 
-Shad charge. From MC simulations it is found to be 
1.008 ± 0.007 for the charged tag and 1.176 ± 0.013 for 
the neutral one. 

The systematic uncertainties associated with the signal 
yield are due to the statistical errors on the PDF param- 
eters (fixed from the MC sample) and potential data/MC 



disagreement for the shapes. We vary the parameters ac- 
cording to their statistical error and correlations. The 
background shapes are validated in the SL (HAD) anal- 
ysis with the mo (ties) sideband: the data/MC ratio of 
the fit variable distribution is parameterized by a first- 
order polynomial, that is used to modify the nominal 
background PDF. A similar strategy is adopted in the SL 
analysis to validate the signal PDF with double-tagged 
events. For the HAD analysis we compare the distribu- 
tions of the NN output before and after the smearing of 
the inputs. In the SL analysis, also the statistical errors 
and the data/MC disagreements for the fractions f s and 
fb are included in these estimates. A summary of the 
systematic uncertainties is listed in Tab. IIVI 

No significant signal is observed in the two analy- 
ses. A Bayesian approach is used to set upper limits 
at the 90% confidence level on B± = B(B+ -> K* + vv), 
Bo = B(B° — > K*°vv) and on their combination. Flat 
prior probabilities are assumed for positive values of both 
B's. Gaussian likelihoods are adopted for the observed 
signal yields, related to the B's by Eq. (P) or Eq. (J5J). 
The Gaussian widths are fixed to the sum in quadrature 
of the statistical and systematic yield errors. We extract 
a posterior two-dimensional PDF P(B±,Bq) using Bayes 
theorem, including in the calculation the effect of sys- 
tematic uncertainties associated with the efficiencies and 
the normalizations, modeled by Gaussian PDFs. Sys- 
tematic uncertainties that are common to the different 
channels and to the two analyses are assumed to be fully 
correlated. The 90% confidence level upper limits are cal- 
culated, after the marginalization of the two-dimensional 
posterior, by: 

/ Po,±(Bo,±)rfB ,± / / Vo,±(Bo,±)dB ,± = 0.9. 

(8) 
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The cross-feed between the different channels is found to 
be negligible in the MC events, but is included in the 
calculation for completeness. We extract the combined 
upper limits: 

B(B+ -> K*+vV) < 8 x 10~ 5 
B(B° -> K*°vV) < 12 x 10~ 5 
B(B -> K*vv) < 8 x 10~ 5 . (9) 

In summary, we search for B — > K*vT> decays in a 
data sample corresponding to 413 fb _1 , collected by the 
BABAR experiment at the T(4S) resonance. We do not 
observe a significant signal in any of the modes studied 
and set upper limits on the decays B° — > K*°vV and 
B + — > K* + vV : and the combined channel B — > K*vv. 
Since no constraints were applied to the kinematics of the 
final state K* meson, or the undetected vv system, these 
results can be interpreted in the context of new physics 
models where invisible particles, other than neutrinos, 
are responsible for the missing energy 0, H, H| . In this 
way, the results presented here are model independent. 
These results represent the most stringent upper limits 
reported to date and they are still consistent with the 
SM expectation [l[. 
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